Instabilities in Binary Mixtures of One-Dimensional Quantum Degenerate Gases 
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We show that one-dimensional binary mixtures of bosons or of a boson and a spin-polarized 
fermion are Luttinger liquids with the following instabilities: i) For different particle densities, 
strong attraction between the mixture components leads to collapse, while strong repulsion leads 
to demixing, ii) For a low-density mixture of two gases of impenetrable bosons (or a spin-polarized 
fermion and an impenetrable boson) of equal densities, the system develops a gap and exhibits en- 
hanced pairing fluctuations when there is attraction between the components. In the boson-fermion 
mixture, the pairing fluctuations occur at finite momentum. Our conclusions apply to mixtures 
both on the continuum and on optical lattices away from integer or fractional commensurability. 
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The exploration of properties of matter at increasingly 
lower temperatures and densities has yielded many sur- 
prises. The achievement of Bose-Einstein condensation 
(BEG) in dilute gases of ultracold atoms (1| is one beau- 
tiful example. Ever since the observation of BEG, one 
line of current experimental (and theoretical) studies is 
on cold gases in very elongated atomic traps. In effec- 
tively one-dimensional (ID) systems strong phase fluc- 
tuations can occur. These systems are interesting for 
both large and small 7 values, where 7 = Mg/h^po is 
the dimensionless coupling ratio, g is the coupling char- 
acterizing the interaction between the atoms in ID, M 
the atom mass, and pa the particle density. For 7^1, 
a single-component ID Bosc fluid is a quasicondensate 
with a fluctuating phase but strongly suppressed den- 
sity fluctuations Large 7 leads to a crossover to 
the Tonks-Giradeau (TG) regime 0: the bosons become 
impenetrable, and resemble noninteracting fermions 0. 
The low-energy properties in both regimes can be well 
described by the Luttinger liquid concept 0,0. 

These aforementioned regimes make up the state space 
of the single-component system of ID repulsive bosons. 
It may therefore seem surprising that adding a second 
component to form a mixture of two types of bosons {B + 
B) or of a boson plus a spin-polarized fermion {B+F) can 
lead to a richer phase diagram. In this Letter, we show 
that such binary mixtures are Luttinger liquids which 
can become unstable to gap-opening in one branch of 
the excitation spectrum, or to demixing, or to collapse. 

These instabilities occur at low temperatures by tun- 
ing some parameters of the system, signalling a quan- 
tum phase transition (QPT) J7]. Some of these instabili- 
ties have direct analogues in higher dimensional systems, 
while others may be more specific to ID systems (and are 
analogues of ID spin- 1/2 fermionic systems [1,0 ). Dilute 
cold gases are ideal for observing the instabilities that we 
predict because, in contrast to conventional solid-state 
materials, they are cleaner and there is a large degree 
of control and tunability of parameters such as particle 
densities and the sign and strength of the interaction. 



Furthermore, diluteness makes it a priori possible to re- 
late microscopic parameters such as scattering lengths 
and atomic masses to the parameters of the low-energy 
effective theory. 

Gontrary to their 3D counterparts 0, UlI IT^ . 
ID binary mixtures have only recently attracted inter- 
ests To our knowledge, there has been no an- 
alytic treatment of these systems within the Luttinger 
liquid framework. This is undertaken here using the har- 
monic fluid approach, which takes into account phase and 
density fluctuations in ID fluids 0,0. We summarize our 
results: 

Case 1.- We consider a binary mixture of two distinct 
atomic species (either B + B ot B + F) or of two internal 
states of the same bosonic species. When a sufficiently 
repulsive interaction exists between the components, the 
mixture will be unstable against demixing, whereas for 
sufficient attraction, it will collapse. 

Case 2.- We consider a binary mixture in a low-density 
limit, when the two components are sufficiently close in 
density and if their sound velocities are similar, the sys- 
tem undergoes a QPT when there is an attraction, no 
matter how weak, between the two components. The re- 
sulting state exhibits a strong pairing tendency. In the 
case oi& B + F mixture, the pairing fluctuations occur at 
the Fermi momentum, unlike the B-\-B case (or F-\-F 0) 
where they occur at zero momentum. 

Model and calculations: The experimental realization 
of the quasi onc-dimcnsional systems requires a tight 
transverse confinement 12, Ijl , which nowadays can be 
achieved in various setups |lftlrRll^ . Our discussion be- 
low will focus on the properties of homogeneous ID sys- 
tems. This should be appropriate for the central region 
of large harmonic 0| or square- well traps 17]. Alterna- 
tively, a toroidal trap 1^ can provide the conditions for 
realizing a homogeneous ID system. When necessary, we 
shall also discuss the effect of longitudinal confinement 
and finite size on the results to be presented below. 

Our starting point will be the following Hamiltonian 
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for a mixture of two ID dilute gases {a,j3 = 1,2): 



H 



\ j dx Pa{x)pi3{x). 



(x) 



(1) 



For bosons, the field operators, '^a{x), obey 
^Jj(a;')] = 6{x — x'), commuting otherwise. 
For a spin-polarized fermion, the field operator obeys 
anti-commuting relations: {5'2(x), ^'3(2;')} — S{x — x'). 
The density operators pa{x) — ^'Jj(x)\E'a(a;). The inter- 
atomic interaction potential is Vaf3{x) = ga/sS^x), with 
gap — 2hLu±aaf} 3], where Oq,^ are the scattering lengths 
parametrizing the 3D interaction between species a and 
(3, and u;± is the transverse confinement frequency. 

To study the low-energy properties of the above model, 
we use the harmonic fluid approach which em- 

ploys the phase-density representation of the field oper- 
ators. Phonon fields are introduced for each component: 
(f>a{x) for long wavelength phase fluctuations and Oa{x) 
for long wavelength density fluctuations. It is assumed 
that dxOa{x) and dx4>a{x) are small compared to the 
equilibrium densities poa- This approach treats bosons 
and fermions in ID on equal footing (^1^, and does not 
assume the existence of a BEC, as usual mean- field treat- 
ments do. Therefore, it can describe both systems of im- 
penetrable bosons (e.g. the Tonks gas Q and others 
and quasicondensates Introducing 
[poa +5:,6'„(x)/7r]^/^X;me™'^'"'°^+^™^°(^)e*'^°(^), where 
the sum over m involves only even (odd) inte- 
gers if the operator is bosonic (fermionic), and 
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[POq + dxOaix)/!!] 



+ 00 



into Eq. Q one obtains the low-energy Hamiltonian: 



7^ {dxdaf + VaKa {dxi'af 
Ka 



%^^^Jdxj2 

— / dx [2gf 8x018x02 + gt cos 2(01 -62 + nSx)] . (2) 
ivr J 



In Eq.(|2J), we have retained only terms which can have 
the most dominant effects at low energies (i.e. marginal 
or relevant in the renormalization group sense). Its va- 
lidity is restricted to energies smaller than the chemi- 
cal potential. The low-energy physics is fully character- 
ized by the phenomenological parameters Ki^2, the sound 
velocities vi^2, the couplings gt^b, and S — poi — P02- 
For small enough gi2, it is possible to relate these phe- 
nomenological parameters to the microscopic ones be- 
cause exchange and correlation effects will be of 0{gf2) 
in gt^b (Otherwise, they must be extracted from nu- 
merics or experiments |y|). First notice that vjya = 
l/{hTrp'^KSa), where Ksa = Pa'^{8pal 8pa) is the com- 
pressibility |5|,|6|. Furthermore, in homogeneous systems, 
Galilean invariance fixes the product vja = v^Ka = 



vpa = hirpQa/Ma Q- Thus from the exact solution for 
a single-component ID Bose fluid VNa = VaK^^ = 
VFa [1 - 87-1 + 0{-f-^)] for 7„ = Magaa/h^Poa > 1 
and VNa — vpalal''^'^ , for 7c < 1. Therefore, 1 < Ka < 
+00 for bosons, Ka = 1 in the Tonks limit. For spin- 
polarized fermions, if 2 ~ 1, because the s-wave scatter- 
ing length vanishes thanks to the Pauli principle, leaving 
the much weaker p-wave channel, which can be neglected. 
Finally, in the weak coupling limit, .g/ — 1712 /(?i7r) and 
gh ~ 47rgi2poiPo2/^- The sign and strength of 1712 can be 
controlled using a Feshbach or a confinement induced Q 
resonance. Interestingly, the same Hamitonian, Eq 13, 
also describes a binary mixture in a ID optical lattice j2l| 
provided that none of the components is commensurate 
with the lattice periodicity. This allows for further pos- 
sibilities of tuning the parameters Va, Ka, and gi2 (e.g. 
the strength of 312 can be modified by shifting relative 
to each other the lattices where each component hops) . 

The above Hamiltonian, Eq. ||2Jl, describes a rich vari- 
ety of one-dimensional binary mixtures: 

Case 1.- When the equilibrium densities of the two 
species are different, i.e. for sufficiently large 5 = poi ~ 
P02, the cosine term in Eq. |5J can be neglected (see be- 
low). Thus one is left with two Luttinger liquids coupled 
by 3/ 8x0x8x92 ■ The normal modes of the system can 
be found from the equations of motion for Oi(x,t) and 
02 (a;, t). The phase velocity of the normal modes is: 

4 - livl + vl) ± y{vl - vlf + A{v.n){vj2)~g}. (3) 
Hence, the mixture will become unstable provided that 



\/VniVN2 < \gf\- 



(4) 



For a repulsive interaction [gf > 0), the two components 
of the mixture repel each other sufficiently strongly to 
demix. However, for an attractive interaction {gf < 0), 
the system is expected to collapse. 

We now examine the linear instability condition, 
Eq. 10}, in several limiting cases. For a, B + B mixture in 
the quasicondensate regime (71^2 ^ 1), VNa — vpaja/'!^'^ 
and gf — gi2l{'hiT). Then the mixture is unstable when 
\/5ii522 < l5i2|- In higher dimensional systems this con- 
dition is obtained from a mean-field theory [l^ . which 
implicitly assumes the existence of a condensate. It 
is therefore interesting that it also applies to ID sys- 
tems in the quasicondensate regime, where the phase 
fluctuates. Another interesting limit is a _B + _F mix- 
ture where the boson is a quasi condensate. For the 
fermion, vn2 — vf2 = TT-T^Pm/^h- Hence, the system 
will become unstable below a critical fermion density 
P02' = -^2'7?9/(7r^? i^qi i). This agrees with the mean- 
field result of Das [lj|. Lastly, consider a B -I- or a 
B + B mixture with impenetrable boson(s) (7^ ^ 1). 
Then Eq. Q implies that the system becomes unstable 
when (P01P02) < MiM2gl2/i'^^^'^)- 
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We now briefly discuss the properties of a stable mix- 
ture. As found above, at long wavelengths there will be 
two branches of modes with phase velocities v± given by 
Eq. ©. This result is not qualitatively affected by the 
presence of a longitudinal confinement, though the actual 
energies of the modes may. More interestingly, a stable 
binary mixture will be a Luttinger liquid with all cor- 
relations characterized (in the thermodynamic limit and 
at zero temperature) by the usual power laws |^. In a 
trapped system, the power laws will be accurate near the 
center of the trap [g. A full discussion of these issues 
will be given elsewhere [2l| . 

Case 2: When the components of the mixture are close 
enough in density, one can no longer neglect the cosine 
term in Eq. To see this, consider the Tonks limit 
where the bosons behave as free fermions. Then this term 
describes a backscattering process between the fermions: 
A fermion of each species is excited from one Fermi point 
to the other, exchanging a momentum close to 2t:{pq±6). 
If 5 <C po = {poi + Po2)/2 then this process can occur at 
low energies and needs to be taken into account. To sim- 
plify our analysis, we take 6 = {S will be addressed 
at the end) . By a perturbative analysis of the free energy, 
we see that the correction due to this term becomes larger 
than the zeroth order term as T — > for Ki + K2 < 2. 
To go further, we use the renormalization group (RG) 
method and take Ki ~ K2 = 1 and vi = V2 (We have 
checked that small deviations from these assumptions do 
not alter the physical picture qualitatively [21|). The 
rationale for this is that the lower bound in a boson sys- 
tem is Ka = 1, which is reached for the Tonks gas limit 
(i.e. by reducing the density, increasing the scattering 
lengths Qaa or both 2]). In a i? + mixture, however, 
the fermion already has K2 — 1, and we only need to 
tune the boson to the Tonks limit. Furthermore, if the 
velocites vi and V2 can be made made similar [22j, then 
in the Tonks limit at equal densities, vi = V2 = vp. Thus 
the Hamiltonian ^ becomes i?off — + iJ_ , where: 



H+ = 
+ 



"2^ 

h 

2^ 



dx 
dx 
dx 



K^\d,e+Y + K+{dM , (5) 



gb cos 



^ix)-gfid^e^f . (6) 



The new fields e± (x) = {9i (x) ± 6*2 (x)) / ^2 and 0± (x) = 
((/>i(a;) ± (/>2(a;)) /\/2 describe the in-phase ("+") and 
out-of-phase ("— ") density and phase fluctuations of the 
components of the mixture. The fact that in this limit 
they are decoupled is an analog of spin-charge separation, 
so ubiquitous in ID Fermi systems |^E3|. 

-ff+ has the form of a Luttinger liquid Hamil- 
tonian characterized by the parameters A'+ = 
[1 + gf /vf]'^^"^ — [1 + .gi2/(?i7ru_F)]~"^^^ and v+ — 
Vp [I + gf /vf]^^'^ — vf [I + gi2/{fi-T^VF)]^^'^ ■ Hence, in- 



phase modes are gapless (in the thermodynamic limit). 

H- is a sine-Gordon model 0. The relative impor- 
tance of the marginal operators, the cosine term and 
{dxO- Y, can be assessed using the RG at weak coupling, 
where 5/ = Qf/vF 1 and gt — p'^'^gb/^i^'^VF ^ 1. 
To second order in these couplings, the scaling flow is 
of Berezinskii-Kosterlitz-Thouless (BKT) type (see e.g. 
Refs. Hi). It is interesting to note that if the sound 
velocities are tuned to be equal (i.e. vi = V2) in the 
limit Ki — K2 = 1, then the RG flow proceeds along 
the separatrix g = gf = gf, on gi2 and the spectrum of 
the system exhibits an enhanced SU(2) symmetry. This 
is not present in the microscopic Hamiltonian and 
seems quite striking, especially for a, B + F mixture. 
Thus, for small gi2 > 0, the effective coupling g is renor- 
malized to zero as the temperature is decreased and the 
cosine term leads only to subleading corrections to the 
low-tcmperature properties: the out-of-phasc modes will 
behave as a Luttinger liquid with = 1 and u_ = vp- 
But when small gi2 < 0, the effective coupling gi, +00 
with decreasing temperature. The exact solution of the 
sine-Gordon model |2| then tells us that a gap opens 
for the excitations of the out-of-phase mode. At weak 
g cx gi2 the gap A ^ l^l^/^e"-^/'^' 0,@1j i-^- exponentially 
small, and it becomes larger as 1712 increases. Hence, at 
low enough temperatures only in-phase modes can be 
excited. In a finite-size system, however, the RG flow to 
strong coupling will be cut-off at a length scale compa- 
rable to the size of the system. Therefore, for weak at- 
traction one should observe that the out-of-phase modes 
become more "stiff" until the full gap develops. 

If vi ^ V2 and/or Ki,K2 ^ 1 (with the differences 
small) the effective SU(2) symmetry is lost and the tran- 
sition is of BKT type, the gap again being exponentially 
small but with a different fimctional form (see e.g. 't','^). 
One notable consequence of the development of the gap 
is that density fluctuations of both components will be- 
come correlated because only in-phase modes are allowed 
at low enough temperatures: specifically, the 2npQ part of 
density correlation {pa{x)p0{O))\ ^ cos(27rpo) x~'^+ 
for arbitrary a — (3 — 1,2 and a (3. These correlations 
can be observed by exciting the modes of the system and 
imaging each species independently. 

The aforementioned correlations are also reflected 
in the "pairing" correlations described by the oper- 
ator A(a;) = "if i{x)'^ 2{x) ■ For a. B + B mixture, 
(A'l'(a::)A(0)) ~ x~^^-^+, which decays more slowly than 
the density correlation, since for gi2 < 0, ~ 

[1 + gi2/{hTrvF)]~^^'^ > 1. For a B -I- F mixture, 
the same correlation function has an interesting oscil- 
lation at wpo (rather than at 27rpo): (^^(2;)A(0)) ~ 
sin {irpox) x~^^^+~^+^^ . Thus, in a B + F mixture, the 
pairing correlations will decay more slowly than the den- 
sity correlations for > 2/^/3 ~ 1.16. 

We now discuss the physical picture underlying the 
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FIG. 1: Schematic phase diagram with the boson(s) in the 
Tonks limit. |poi — P02i is the density difference at (/12 = 0. 



above calculations. On intuitive grounds, an attraction 
between the two components of the mixture {g\2 < 0) will 
lead to the formation of bound pairs, and might eventu- 
ally lead to collapse. However, in the regime studied 
above, pairing takes place between hard-core bosons or 
a hard-core boson and a fermion, which means that as 
long as the repulsion between bosons of the same species 
is stronger, only bound pairs containing one particle of 
each component are possible. This picture is confirmed 
by considering a strong-coupling limit of the Hamilto- 
nian where first we take gii,(?22 +00 and then 
512 —00 such that \gvi\ ^ 511,522 H^l- Moreover, in 
a B + B mixture pairs are also bosons, but their con- 
densation is forbidden in ID dimension by rigorous the- 
orems 2j. In a _B -|- _F mixture, the pairs are fermions, 
and this slightly suppresses pairing correlations relative 
to the B + B case. We expect a smooth crossover from 
the weak coupling regime to the strong-coupling regime 
where tightly bound pairs form. As the strength of at- 
traction between the two components increases, the size 
of the gap also grows thus preventing collapse (Fig.QJ. 

However, the gap can be destroyed when a sufficiently 
large difference in chemical potential leads to a density 
difference \S\ > 0. This transition is of commensurate- 
incommensurate type 0. Small inhomogeneities in the 
density due to the trap, will not have a strong effect in the 
pairing provided that the density difference between the 
two components does not become too large. Otherwise, 
one can expect phase segregation into pairing and non- 
pairing regions, with the former occurring mainly near 
the center of the trap. 

In conclusion, we have shown that binary mixtures of 
ID quantum gases exhibit some interesting instabilities. 
In particular, we have found a pairing instability which 
may be special to one dimension and should have unusual 
experimental signatures [2ll |. 

We acknowledge inspiring conversations with A. Ners- 
esyan and M. Fabrizio, and useful correspondence with B. 
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